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Stability of continuously pumped atom lasers 
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S. A. Haine, J. J. Hope, N. P. Robins and C. M. Savage 
Department of Physics and Theoretical Physics, Australian National University, ACT 0200, Australi^ 

A multimode model of a continuously pumped atom laser is shown to be unstable below a critical 
value of the scattering length. Above the critical scattering length, the atom laser reaches a steady 
state, the stability of which increases with pumping. Below this limit the laser does not reach a 
steady state. This instability results from the competition between gain and loss for the excited 
states of the lasing mode. It will determine a fundamental limit for the linewidth of an atom laser 
beam. 

PACS numbers: PACS numbers: 03.75.Fi, 03.75.Be 



Introduction. -The physics of lasers, both fundamental 
and applied, has been an important theme in modern 
physics. Recently, the development of dilute gas Bose- 
Einstein condensation jjj has opened up the study of 
atom lasers: the matter wave analogues of optical lasers 
0). Many of the defining characteristics of optical lasers, 
such as high spectral density, recur for atom lasers. In 
this letter we report on a critical, and somewhat counter- 
intuitive, difference between atom and optical lasers. One 
might expect that an atom laser would behave most like 
an ideal optical laser when the two-body interactions be- 
tween atoms were minimal, because this is true for pho- 
tons in an optical laser. We show that this is not the case, 
and that stable atom laser operation requires sufficiently 
repulsive two-body interactions between the atoms. That 
these interactions are necessary is particularly significant 
because they induce a phase diffusion in the lasing mode, 
which places a fundamental limit on the linewidth of the 
atom laser beam. 

Outcoupling from magnetically trapped BEC has been 
achieved by using radio frequency radiation |], |J] or Ra- 
man transitions H to change the internal state of the 
atoms, causing them to be either untrapped or anti- 
trapped. Interference has been observed in the output 
field H , and quasi-continuous beams have been produced 
j| ||. Coupling the atoms out more slowly causes the 
output beam to have a narrower linewidth, but this comes 
at the expense of the beam flux |Q. High spectral flux 
in optical lasers is generated through a competition be- 
tween a depletable pumping mechanism and the damping 
of the lasing mode. This means that a higher pumping 
rate leads not only to greater flux, but also to a narrower 
output spectrum Q . An atom laser with gain- narrowing 
must also have a saturable pumping mechanism that op- 
erates at the same time as the damping B. 

Atom lasers have been described semiclassically with 
a large number of modes, and with a full quantum me- 
chanical model containing a few modes. Multimode mod- 
els, using the Gross-Pitaevskii equation M, are capable 
of including the spatial effects of the atomic field and 
the interatomic interactions p!o| . This makes them use- 
ful for describing non-pumped atom laser experiments in 



which these effects are dominant. Multimode models are 
required to determine whether pumped lasers approach 
single mode operation. As for optical lasers, if the lasing 
system can be accurately described as a single mode then 
the quantum statistics of that mode control the linewidth 
of the output beam. This means that although a semi- 
classical model is useful to determine the stability of an 
atom laser, a full quantum mechanical model must be 
used to calculate the linewidth of the beam coming from 
a pumped atom laser. Single (or few) mode quantum 
theories have been used to describe continuously pumped 
atom lasers O] . Unfortunately, these models cannot de- 
scribe the multimode behaviour of the atom laser such 
as the non-Mar kovian nature of the output coupling [jL2| , 
except under the approximation that the lasing mode 
itself is single mode [(L3). It is also extremely difficult 
to include the effects of the atomic interactions in these 
models. Atom laser linewidths will be limited by these in- 
teractions due to Kerr-like dephasing of the lasing mode 
|l4| , as well as by thermal effects ]lf| . 

Including the spatial effects, interactions and quan- 
tum statistics of the atom field requires a quantum field 
theory calculation which may become tractable by using 
stochastic field theories derived from phase space meth- 
ods [[l6) . These stochastic field simulations reduce to the 
semiclassical model in the limit of zero noise, so the semi- 
classical stability of the relevant system must first be ex- 
amined. We describe a semiclassical, multimode model of 
an atom laser and investigate the effects of varying the 
strength of the atom-atom interactions relative to the 
trap strength. This parameter can be controlled experi- 
mentally by using Feshbach resonances |l7j or by adjust- 
ing the trap parameters. We find that the atom laser is 
stable only when the repulsive interactions between the 
lasing mode atoms are sufficiently strong. For a given 
pumping rate and trap parameters, this results in the 
model being stable above a critical scattering length and 
unstable otherwise. 

Atom laser modeL-Although a full quantum field the- 
ory will be required to calculate the linewidth of an atom 
laser with interactions, the stability of such a device can 
be determined with a multimode semiclassical treatment. 
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We model the atom laser as a two-component semiclas- 
sical atomic field coupled to an incoherent reservoir of 
atoms described by a density p(x). The first field, ipt( x )> 
is trapped in a harmonic potential and will form the las- 
ing mode. The second field, ip u (*-), is untrapped and 
will form the laser beam. The coupling can be achieved 
either by radio frequency transitions between the elec- 
tronic states H or by a Raman transition ||. The 
fields obey equations of the Gross-Pitaevskii form, with 
damping and reservoir coupling. The dynamical equa- 
tions for the coupled system are 

■dipt , h 2 V . (1) U tu . |2 , U u , , ,2 

-£ = r - JpP- K p \ipt\ 2 p + AV 2 p, (1) 

where M is the mass of the atom, V is the trapping 
potential, g is the acceleration due to gravity, which is 
assumed to be directed along the x direction. = 
4irha,ij/M is the interatomic interaction between ipi and 
tpj , dij is the s-wave scattering length between those same 
fields, 7j is the loss rate of ipi due to background gas 
collisions, 7 4 is the loss rate of ipi due to two-body in- 
elastic collisions between particles in that state, j tu ' is 
the loss rate of each field due to two-body inelastic colli- 
sions between particles in the other electronic state, kr 
is the coupling rate between the trapped field and the 
output beam, k is the momentum kick due to the cou- 
pling process, k p is the coupling rate between the pump 
reservoir and the trapped field. For the pump reservoir, 
r is the rate of density increase, 7 P is the loss rate, and 
A is the spatial diffusion coefficient. 

We have not included three-body losses, which are nor- 
mally important. However, near a Feshbach resonance 
they may be negligible, as for 85 Rb, and still allow a 
wide range for the scattering length |l^] . Our results are 
not affected by this limitation, but since it is not known 
how to include three-body losses in stochastic field simu- 
lations, we ignore them to facilitate later extension to a 
quantum field model. 

The pumping terms in the above equations are phe- 
nomenological, describing an irreversible pumping mech- 
anism from a reservoir which can be depleted, but is re- 
plenished at a steady rate. These two features are nec- 
essary for any pumping mechanism that generates gain- 
narrowing through the competition of the gain and loss 
processes of the lasing mode. 

It has been shown that output coupling can induce 
a localised bound state jl^, [l3|. When the atom laser 
is pumped, the population of this bound state increases 



indefinitely. Gravity must be included in the model to 
make this state decay, and allow steady operation to be 
achieved. For a given set of parameters, there is a maxi- 
mum coupling rate kr below which this metastable state 
has negligible effect on the total dynamics. We have also 
included the inelastic losses, which improve the overall 
stability. 

Stability of the atom laser modeZ.-Throughout the rest 
of this letter we analyse the model given in Eq.(|l|) in one 
dimension, assuming that both the lasing mode and the 
output beam are strongly trapped in the transverse di- 
mensions. We use the atomic properties and loss rates for 
85 Rb near a Feshbach resonance, where the interatomic 
interactions can be tuned with the magnetic fields [jl7] . 
We will use the following parameters for all subsequent 
calculations: 7t (1) = 7.0 x ICT 3 s~\ 7t (2) = 1.7 x ICT 8 
ms- 1 , 7 i 2) = 3.3 x ICT 9 ins" 1 , 7t (2) = 0.5 7t (2) , 7p = 5 
s~ x , k p — 6.3 x 10~ 4 ms -1 , a trapping potential of 
V = mu} 2 x 2 /2, where w = 50 rad s _1 , |k| = 10 6 m _1 , 
A = 0.01 m 2 s _1 . The interatomic interactions between 
different Zeeman levels are unknown for 85 Rb, so we as- 
sume: U t = U u = 2U tu , a choice which has negligible 
effect on the results of this work. 

We integrate Eqs.(Q) numerically using a psuedo- 
spectral method with a Runge-Kutta time step It 
was found that the atom laser model exhibits varying 
levels of stability depending on both pumping rate and 
scattering length. Figure |l| shows the density of the atom 
laser beam at a point below the condensate for various 
pumping rates, and for two scattering lengths. The total 
pumping rate is R = rL, where L is the width of the 
pumping region. We find that for zero scattering length 
the laser is unstable, with increasing pumping leading 
to modal oscillations. In contrast, for large scattering 
length the atom laser becomes stable above a threshold 
pumping rate, and increasingly stable as the pumping 
rate increases. 

The dominant oscillation of the laser mode is a breath- 
ing mode. As we shall see, increasing the pumping in- 
creases the damping of this, and the higher frequency 
modes, but only at a high scattering length. Varying the 
scattering length between the two cases shown in Fig. 1 
gives intermediate types of behaviour, ranging from un- 
stable to stable, with an increasing scattering length in 
the stable case making the oscillations more strongly 
damped. 

Analysis of the a u — case. -With two-body interac- 
tions off, the spatial shape of the lasing mode is indepen- 
dent of the number of atoms in the trap. An atom laser 
operating in this regime experiences no nonlinear phase 
diffusion, but then our atom laser model is not stable. 
The semiclassical instability of this system can be under- 
stood with a simplified model, as it is directly related to 
the eigenfunctions of the lasing mode. In this model we 
ignore the output coupling and hence the equation for 
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FIG. 1: Density of the atom laser beam (in mm -1 ), at a point 
below the trapped region, as a function of time, (a) Scattering 
length att = and coupling rate kr = 10 s -1 . Pumping rates 
are (i) R = 10 5 s _1 and (ii) R = 10 6 s" 1 . (b) a tt = 20 nm 
and kr = 300 s -1 . Pumping rates are (i) 7? = 500s -1 , (ii) 
R = 10 4 s -1 , (iii) R = 10 5 s -1 , and (iv) R = 10 6 s" 1 . All 
other parameters are as given in the text. 



the laser mode becomes 



. dipt 
dt 



T + V 



i7t (1) - i7t (2 V t | 2 + ^ P AyV> t , (2) 



where T is the kinetic energy operator, and N p — 
pdx. Expanding Eq.(^) into the basis of harmonic 

oscillator eigenfunctions {<fi n (x) : (T + V)<p n = E n <j) n } 
and solving for the mode coefficients gives the equivalent 
set of equations 
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with I mjk = J™ oa (l>n<i>i<t>j<t>kdx and E n = (n + 
and where we interpret |c„| 2 as the number of atoms in 
mode n. The two-body loss term serves as a form of cou- 
pling between modes of equal parity, with the diagonal 
terms /„ = I n nnn acting as damping for each mode. The 
first few diagonal terms are Iq = 0.5/3, Ji = 0.375/3, 12 = 
0.320/3, h = 0.287/3, h = 0.264/3, 1 5 = 0.244/3 with 
/3 = VA^FH^)" They decrease with increasing n, so the 
excited states experience lower loss than the lower energy 
states. The loss rate decreases because the excited states 
are more spread out than the lower energy states. This 
same trend occurs for three-body loss. The pumping is 
not mode selective, so the excited states have a higher net 
gain. The coupling between modes seeds excited states 
which become more populated than lower states. This 
growth is the origin of the atom laser instability. 

Stability of a pumped Bose-Einstein condensate.-With 
two-body interactions on, the previous simplified model 
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FIG. 2: Time series and power spectra of the laser mode 
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scattering lengths are (i) a t t = 0.01 nm, (ii) a t t = 0.5 nm, 
(iii) a t t = 1 nm, and (iv) att = 10 nm. Spectral power is in 
arbitrary units and is obtained from the time series after the 
initial 0.4 s. 



is inadequate. As the interactions between the atoms be- 
come dominant, the excited states have a different spatial 
shape, and may become more damped than the ground 
state. For mode-insensitive pumping this would lead to 
a selection of the ground state. 

Hence we expect there to be a critical scattering length 
for which the atom laser moves from unstable to stable. 
Since the nonlinear interactions depend on the density, 
we also expect the critical scattering length to be depen- 
dent on the pump rate and trapping frequency. 

We found that the semiclassical dynamics of the atom 
laser beam followed the dynamics of the lasing mode con- 
densate. For the range of output coupling rates that we 
consider, the condensate dynamics is only weakly affected 
by the output coupled beam. Hence, in order to reduce 
the computational load, we ignored the output coupling 
in the following. 

Figure || shows the different regimes of stability as the 
scattering length is changed with fixed pumping rate. 
The case with the lowest scattering length is unstable 
due to the excitation of higher trap modes. The sys- 
tem with the highest scattering length is stable, with all 
excited states being damped. The systems with inter- 
mediate scattering length show a reasonably steady pop- 
ulation in the lowest excited states, and a slow growth 
of high frequencies due to the population of the excited 
states. Figure ^ also shows the power spectrum of the 
oscillations in the central density, ignoring the transients 
in the first 0.4 seconds of each trace. The spectrum of the 
case with the lowest scattering length shows that the be- 
haviour resembles that of the harmonic oscillator, with 
growing oscillations at the frequencies 2nuj = 100n/2vT 
Hz, which are the even harmonic oscillator eigenfrequen- 
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FIG. 3: Atom laser stability in pumping rate and scattering 
length parameter space (log- log scales). Cases which we have 
determined to be in the stable region (Type I) are shown as 
squares, Type II systems are shown as circles, and unstable 
(Type III) systems are shown as stars. The bold lines separate 
the stability regions. 



cies. At larger scattering lengths, the behaviour is stable, 
but the transient oscillations closely resemble the even 
Thomas- Fermi eigenfrequencies uj^Jn(n + l)/2 (n even) 

We found it convenient to classify the behaviour into 
three catogories. Type I behaviour is where all excited 
modes are damped. In Type II behaviour the breath- 
ing mode is damped but higher modes are undamped, 
so oscillations are initially decreasing, but the system is 
unstable over a long period. In Type III behaviour all 
modes are undamped, and the system is unstable. Fig- 
ure ^ shows the type of behaviour found for different 
values of the pumping rate R and scattering length a«. 
For high pumping rates, the high frequency oscillations 
produce numerical instabilities, so the boundaries in the 
phase diagram were difficult to explore in detail. 

The semiclassical stability of the pumped condensate 
is almost purely determined by the pumping rate and 
the scattering length, as these two parameters control 
the mean field interaction strength. 

Conclusions: We have shown that a semiclassical 
atom laser becomes stable above a critical scattering 
length. This is because the excited states are increas- 
ingly damped when the interactions are the dominant 
contribution to the lasing mode energy, and otherwise 



they are more weakly damped than the lower energy las- 
ing modes. This effect is not dependent on the output 
coupling, but could possibly be controlled by a highly 
mode selective pumping mechanism. This result is im- 
portant in the design of continuously pumped atom lasers 
because the Kerr-like phase diffusion will scale with the 
interactions, and will induce a fundamental limit to the 
linewidth of the atom laser beam |l4j . Due to the oppo- 
site scaling of the semiclassical and the quantum noise, if 
an atom laser can be operated in a stable regime then we 
expect it to have a minimum linewidth at an optimum 
scattering length. 

This research was supported by the Australian Re- 
search Council and by an award under the Merit Allo- 
cation Scheme on the National Facility of the Australian 
Partnership for Advanced Computing. 



[3] 

[5] 
[6] 
[7] 
[8] 

[9] 
[10] 

[11] 



[12] 
[13] 
[14] 

[15] 
[16] 
[17] 
[18] 



Electronic address: joseph.hope@anu.edu.au] 

F. Dalfovo et al, Rev. Mod. Phys. 71, 463 (1999); A. J. 
Leggett, Rev. Mod. Phys. 73, 307 (2001). 

R. J. Ballagh and C. M. Savage in Bose-Einstein Con- 
densation, edited by C. M. Savage and M. P. Das (World 
Scientific, Singapore, 2000). 

I. Bloch et ai, Phys. Rev. Lett. 82, 3008 (1999). 

J. L. Martin et al, J. Phys. B 32, 3065 (1999); M.-O. 

Mewes et ai, Phys. Rev. Lett. 78, 582 (1997). 

E. W. Hagley et ai, Science 283, 1706 (1997). 

I. Bloch et al, Nature 403, 166 (2000). 

J. J. Hope, Phys. Rev. A 55, 2531 (1997). 

L. Mandel and E. Wolf. Optical Coherence and Quantum 

Optics (Cambridge University Press, Cambridge, 1995). 

H. M. Wiseman, Phys. Rev. A 56, 2068 (1997). 

N. Robins et al, Phys. Rev. A 64, 043605 (2001). 

H. M. Wiseman and M. J. Collett, Phys. Lett. A 202, 

246 (1995); H. Wiseman et ai, Quant. Semiclass. Opt. 8, 

737 (1996); M. Holland et ai, Phys. Rev. A 54, R1757 

(1996). 

G. M. Moy et al, Phys. Rev. A 59, 667 (1999). 

J. J. Hope et ai, Phys. Rev. A 61, 023603 (2000). 

H. M. Wiseman and L. K. Thomsen, Phys. Rev. Lett. 
86, 1143 (2001). 

R. Graham, Phys. Rev. Lett. 81, 5262 (1998). 
M. J. Steel et ai, Phys. Rev. A 58, 4824 (1998). 
J. L. Roberts et ai, Phys. Rev. Lett. 85, 728 (2000). 
The method, known as RK4IP, has been developed 
by the BEC theory group of R. Ballagh at the Uni- 
versity of Otago. It is described in the Ph.D. the- 
sis of B. M. Caradoc-Davies which is online at: 



http://www.physics.otago.ac.nz/bec2/bmcd/ 



[19] B. Kneer et ai, Phys. Rev. A 58, 4841 (1998). 



